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We report an analysis of neutron diffraction from single crystals of the spin-liquid pyrochlore 
Tb2Ti207 under the application of magnetic fields along the crystallographic [110] direction. Such 
a perturbation has been shown to destroy the spin liquid ground state and induce long-range order, 
although the nature of the ordered state was not immediately determined. Recently, it has been 
proposed that the ordered state is characterized by spin-ice-like correlations, evincing an emergent 
ferromagnetic tendency in this material despite the large negative Curie-Weiss constant. Here, 
we argue instead that the ordered state is dominated by Q / correlations that emerge either 
from strong antiferromagnetism or magnetoelastic distortion of the crystal. In contrast to previous 
reports, we observe no evidence for re-entrant behaviour in the high field limit. Extreme sensitivity of 
the ordered state to the alignment of the applied field is suggested to account for these discrepancies. 



PACS numbers: 75.25.+Z, 75.30.Gw, 75.30.Kz 

As originally conceived by Anderson^ a spin liquid 
is a dynamically disordered quantum antiferromagnetic 
(AFM) phase wherein strong magnetic correlations ex- 
ist, but do not extend beyond near-neighbour distances. 
This behaviour is dependent on geometrical frustration 
in the crystal lattice, where the local connectivity of the 
magnetic ions prohibits the simultaneous satisfaction of 
all near-neighbour interactions? 2 - This frustration is com- 
monly realized by a combination of triangular or tetra- 
hedral architectures and AFM interactions. Ferromag- 
netic interactions can also be frustrated in the presence 
of strong single- ion anisotropy, as occurs in the spin ices? 3 - 

The search for materials that exhibit the ideal spin 
liquid physics envisioned by Anderson has occupied ex- 
perimentalists for decades, with mixed results? 4 - Typi- 
cally, these efforts have focussed on the highly frustrated 
kagome and pyrochlore lattices, in two and three dimen- 
sions respectively. In general, it has been difficult to 
identify a material with spin 1/2 moments that crystal- 
lizes into a perfect kagome or perfect pyrochlore lattice, 
where only near- neighbour couplings are relevant. How- 
ever, a host of interesting materials with slightly more 
complicated interactions have been identified, and they 
exhibit physics that is intriguing in its own right? 2 - 

Tb2Ti2 07 belongs to the family of rare-earth- metal 
titanate pyrochlores, which have garnered recent atten- 
tion due to the exotic frustration-driven physics they 
manifest! 5 - Tb2Ti207 is isostructural to the canonical 
spin ices Ho2Ti2C-7 and Dy2Ti207,— with magnetic rare- 
earth ions occupying the pyrochlore lattice (Fig. 1). 
In contrast to these frustrated ferromagnets, Tb2Ti207 
presents a significantly reduced single-ion anisotropy 
(A ~ 18 K) as well as a negative Curie- Weiss constant 
of 0cw ~ — 19 K? 6 - Interest in Tb2Ti207 was ignited 
when it was observed that dynamic short-range mag- 
netic correlations appear on cooling through T ~ Ocwi 



but that long range magnetic order is absent down to the 
lowest measurable temperatures (~ 20 mK)X These are 
the key experimental signatures of a spin liquid ground 
state. However, Tb2Ti 2 C>7 is expected to be governed 
by a more complicated Hamiltonian than those gener- 
ally considered to model spin liquids, since the magnetic 
moments of the Tb 3+ ions are large and strongly influ- 
enced by spin-orbit coupling and the crystal electric field. 
One would therefore expect the moments in Tb2Ti207 
to be reduced from the free ion value, and to be ori- 
ented along local <111> directions that join the centers 
of the corner-sharing tetrahedra in the latticed This has 
been experimentally confirmed?^ and yet a puzzle re- 
mains: naively, one would expect single- ion anisotropy 
to lift the degeneracy of the spin liquid ground state, and 
drive the system to a non-colinear ordered Neel state? 6 -*! - 
The manner in which frustration is reintroduced is an 
open topic, although recent work has shown that quan- 
tum fluctuations out of the Ising-like ground state dou- 
blet can renormalize the effective interactions to be fer- 
romagnetic, leading to a fluctuating quantum spin ice 
state at low temperatures?! 1 - There also exists experi- 
mental evidence pointing to magnetoelastic effects at low 
temperatures ; 12 ' 13 which suggests that spin-lattice cou- 
pling is an important ingredient of the ground state. 

Regardless of the exact origin of the exotic ground 
state in Tb2Ti207, it is now well known experimen- 
tally that long range order can be induced by the ap- 
plication of magnetic fields^!^ or pressure?! 5 - It is not 
unreasonable to expect that an understanding of these 
perturbation-induced ordered phases could unravel some 
of the mysteries of the unperturbed ground state. There- 
fore, a careful characterization of the ordered states is 
crucial. The original discovery of magnetic field induced 
order— did not collect enough magnetic Bragg reflec- 
tions to facilitate a full magnetic structure refinement, 
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FIG. 1: (Colour online) The pyrochlore lattice. The four 
tetrahedra in the unit cell are labeled by (A,B,C,D) and the 
four magnetic sublattices are labeled by (1,2,3,4). Chains 
running along the [110] applied magnetic field direction are 
known as a-chains, and are shown in red. Perpendicular j3- 
chains are shown in blue. A classical expectation for AFM 
Ising spins in high field is shown. 



TABLE I: Phase factors in the 16 spin unit cell 
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although it did identify the high-field ordering wavevec- 
tor as Q — (27r/a)(112), where a sharp, intense Bragg 
reflection was observed at a Q ^ position. A weaker, 
diffuse feature centered at (112) was also observed, which 
forms a halo around the sharp peak. This is evident in 
Fig. lc in Ref. [13. Subsequent studies involving many 
more magnetic reflections returned qualitatively differ- 
ent results^ finding instead that the dominant scatter- 
ing was at the Q — positions. This was interpreted 
as evidence for the predicted spin ice like character of 
Tb2Ti2C>7. These measurements found weaker evidence 



for the Q ^ phase, with a reduced volume fraction 
and a transition temperature that was suppressed at 
high fields 8 relative to the original work. This phase 
was refined as a polarization of spins with Ising axes 
that project along the magnetic field (generally referred 
to as a-chains), and AFM ordering of spins with Ising 
axes perpendicular to the magnetic field (known as (3- 
chains). A depiction of this phase with a full volume 
fraction, which is an AFM analogue to the Q = X phase 
in the spin iceS ) 16 i 17 is illustrated in Fig. 1. Assum- 
ing comparable sample quality, the main difference be- 
tween the original measurements of Rule et ali^ and 
the subsequent measurements of Cao et o/4 is that the 
original experiment was performed with the magnetic 
field aligned along [110] to within an accuracy of ±0.5°, 
within the mosaic spread of the single crystal sample, 
while the subsequent experiment was performed with the 
magnetic field misset from [110] by ~ 5° according to the 
authors^ A magnetic field misalignment of this order 
has been shown to qualitatively change the ground state 
in the sister spin-ice compounds, both experimentally 18 
and theoretically— Therefore, it is interesting to investi- 
gate the discrepancies between these two measurements 
in further detail. 

Here we report a careful analysis of Bragg reflec- 
tions collected using time-of-flight neutron scattering, in 
the case where the magnetic field is nominally perfectly 
aligned (i.e. to within the mosaic spread of the crystal). 
We present parametric studies of the nuclear allowed 
(113) peak and the superlattice (002) and (112) peaks. 
Although insufficient for a brute-force magnetic struc- 
ture refinement, the relative intensities of these reflec- 
tions strongly constrain the nature of the ordered state. 
Measurements were performed using the Disk Chopper 
Spectrometer at the NIST Center for Neutron Research, 
with cold (5A) incident neutrons^ This time-of-flight 
instrument is immune to higher-order wavelength con- 
tamination problems systematic to diffractomctcrs cm- 
ploying a crystal monochromator. Bragg intensities were 
collected as a function of temperature and magnetic field 
using an 11.5 Tesla superconducting magnet system with 
a dilution insert, and normalized to incoherent vanadium 
scattering to correct for relative detector efficiencies and 
cryostat dark angles. The results are shown in Fig. 2, 
where we plot the results of the raw measurement in 
2a and 2b, and a useful rescaling of the data in 2c and 
2d. This rescaling is performed as follows: the zero- 
field component is subtracted out, and intensities are 
scaled to incoherent vanadium scattering to correct for 
efficiency. Then, we multiply by the magnitude of the 
momentum transfer, take the square root, and divide by 
the magnetic form factor calculated in the dipole approx- 
imation. The resulting quantity is directly proportional 
to the total scattering length in the unit cell which gives 
rise to the Bragg peak, and therefore also proportional 
to the related ordered magnetic moment transverse to 
Q. It is clear from both the raw and rescaled intensities 
that the dominant scattering in the field-induced ordered 
phase occurs at the superlattice (112) position. It is im- 
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TABLE II: Ordered moment in the unit cell. (0.4K, 7T). 
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FIG. 2: (Colour online) Temperature and applied magnetic 
field dependence of Bragg scattering at three reciprocal space 
positions. 2a and 2b show integrated intensities measured by 
rocking the crystal through the Bragg angle, as a function of 
magnetic field and temperature. 2c and 2d show a rescaling 
of this data as described in the text. The low-field hysteresis 
shown in (c) has been previously discussed.— The remnant 
intensity at high temperature at the (ll2) position is due to 
the diffuse halo discussed in the text. 



portant to note that the intensity of this Bragg peak 
continues to increase up to the highest applied magnetic 
field (9 Tesla), in contrast to the re-entrant behaviour 
recently reported to set in above ~ 5 Tcsla4 The two 
preceding points illustrate a clear and qualitative dif- 
ference between the results presented here, and those 
reported by Cao et o/4 They report weak scattering at 
the (112) position that is further suppressed for large 
applied magnetic fields; we report exactly the opposite. 

It is possible to assign qualitative significance to the 
intensities of each of the measured Bragg reflections, if 
we assume that the scattering is magnetic in origin. Ta- 
ble 1 shows the phase factor e l ® r calculated for each 
of the three reflections, and for each spin in the cubic 
nuclear unit cell. The magnetic ion sites are labeled to 
coincide with Fig. 1, where each of the four unit cell 
tetrahedra (A,B,C,D) contain four spins with different 
easy-axes (1,2,3,4). It becomes immediately clear that 
ferromagnetic ordering or field-induced polarization will 
not contribute to the scattering at (112), and that AFM 
correlations between (A,D) tetrahedra and (B,C) tetra- 
hedra will. Further, an AFM moment on the /3-chains 
does not interfere with an AFM moment on the a-chains 
- the two contributions are decoupled since the former 
contributes to the imaginary part of the cross-section, 
while the latter contributes to the real part. It is also 
evident that a net polarization of a-chains will contribute 
to scattering at (113) but a net polarization of /3-chains 
will not, and that (002) scattering is a measure of the 
vector difference in net polarization between the a- and 
/3-chains. Therefore, despite the fact that this subset 
of reflections is insufficient to fully refine the magnetic 
structure, we can confidently comment on the size of the 
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AFM ordered moment, the difference in polarization be- 
tween the two sets of chains, and the degree to which the 
a-chains are polarized along the applied magnetic field. 
By normalizing to the nuclear scattering at (113) mea- 
sured and high and low temperature in zero field, we can 
extract a lower limit on the sizes of these moments which 
is constrained only by an incomplete understanding of 
extinction effects in our crystal, which we neccessarily 
ignore. 

Let us begin by considering the AFM analogue of the 
Q = X spin ice phasei&l 7 . (QXAFM phase), which is 
shown in Fig. 1. This is a classical ground state for 
an Ising AFM in high field, where field perpendicular 
spins (on /3-chains) order into AFM chains, with a phase 
difference of tt between the (A,D) tetrahedra and the 
(B,C) tetrahedra, while a-chain spins maximize their 
polarization along the field. Although it is possible to 
construct a, Q — state in this way as well, where (A,D) 
/3-chain spins are not flipped relative to (B,C) /3-chain 
spins, such a configuration would not generate scatter- 
ing at (112). The classical expectation would be that the 
states described above are degenerate for nearest negh- 
bour interactions, so the selection of Q ^ order is in- 
dicative of the importance of long-range interactions. A 
ground state resembling QXAFM was refined as a minor- 
ity phase in the work of Cao et al.^ but here we consider 
a full volume fraction. By considering the phase factors 
displayed in Table 1, it is clear that such a ground state 
would have identical magnetic scattering cross-sections 
at (002) and (113), which would be proportional to the 
net polarization of the a-chain spins. An inspection of 
Fig. 2 will satisfy that this is not the case experimen- 
tally. The reduction of (002) with respect to (113) im- 
plies some polarization of the beta chains along the field, 
away from their easy-axes. The cross-section at (112) is 
due solely to the component of the /3-chain spins which 
is AFM ordered and transverse to Q, which is reduced 
from the single- ion /3-chain moment by a factor of 2/3. 
Therefore, even by promoting the minority phase of Cao 
et al. to a full volume fraction, we do not recover enough 
scattering power at (112) to account for the measured in- 
tensity, assuming that a and /3 chain moments have the 
same magnitude in the ground state. One would expect 
the scattering at (002) and (113) to be dominant in this 
case. Worse still, we have the unfortunate problem of 
requiring a larger than expected AFM moment on the 
/3-chains while at the same time requiring a finite polar- 
ization of these same spins. Clearly, the classical Ising 
expectation is not a good approximation to the field- 
induced ground state realized experimentally. 

We can better characterize the nature of the ordered 
state by relaxing the Ising constraint. The necessity 
of this procedure is not unexpected, since the rela- 
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tively small anisotropy gap is expected to allow for some 
restoration of rotational symmetry ^ 10 ' 11 In general, the 
symmetry of the problem suggests that we consider in- 
dependent [110]-field-polarized moments on the a and 
(3 chain sites, which we will refer to as P a and Pp re- 
spectively. In addition, we consider independent AFM 
ordered moments on each sets of chains: S a and Sp, 
which contribute to the scattering at (112). In this pic- 
ture, the scattering at (002) will arise from \8(P a — Pp)\, 
the scattering at (113) will arise from |8P Q |, and the scat- 
tering at (112) will probe the quantity \8(S a + iSp)\±. It 
is important to note that while the polarized moments 
are transverse to Q since the magnetic field is applied 
transverse to the scattering plane, for the AFM moments 
this is not necessarily the case. Only the AFM moment 
transverse to Q will contribute to the scattering at (112). 
In order to convert the quantities plotted in Fig. 2c, d 
to Bohr magnetons, we first normalize to the nuclear 
contribution measured at (113), for which the scattering 
length can be calculated for the known crystal lattice in 
cm~ 12 . This value is then converted to Bohr magnetons 
using the ratio The resulting total moments 

transverse to Q, summed over the 16 spin unit cell, are 
listed in Table 2. These values constitute lower limits, as 
extinction effects will act to decrease the measured mag- 
netic moments. However, since extinction is expected 
to preferentially suppress strong reflections, one would 
expect that the dominant effect of a proper extinction 
correction would be to increase the AFM ordered mo- 
ment relative to the polarization. 

Now, we are in a position to offer a description of the 
ordered state. The net polarization per spin, given by 
the quantity (P a + Pp)/2, can be derived as ~ 1.7 if 
P a > P/3, and ~ 3.5 (ib if Pa < Pp- These values are 
significantly lower than the saturation magnetization re- 
ported in the literature 4»22. which could either indicate 
that a strong extinction correction is necessary, that pre- 
cise alignment of the applied magnetic field leads to a 
reduction in the polarized moment, or both. As for the 
AFM ordered moment, if we assume that it is entirely 
transverse to Q, we arrive at an average value of 3.7 /ib 
per spin. However, this represents an extreme departure 
from Ising anisotropy. If we instead assume that the 
AFM moment remains Ising-like, we arrive at a value 
of of 4.9 fiB per spin. Implicitly, we have assumed that 



the moment is equally distributed amongst all spins (i.e. 
| S a | = | Sp | = S). This is neither necessary or likely. 
However, we do so because it is an instructive way to 
estimate the average value. If we assume that only the 
/3-chain spins participate in the AFM order, then we find 
that \Sp\ = 7.8 un- 
it bears mentioning that the above discussion of the 
ordered state, like that of Cao et a/4, is based on unpo- 
larized neutron scattering measurements. It is therefore 
possible that the nuclear Bragg scattering could man- 
ifest some magnetic field dependence, which would be 
impossible to deconvolute from the response of the mag- 
netic spins in these measurements. Tb2Ti20y is a known 
giant magnetostrictive material^, although no evidence 
for a magnetic field induced structural phase transition 
has ever been reported. Nevertheless magnetoelastic dis- 
tortion of the lattice has not been ruled out, and may 
well play a role in the nature of the ordered phase. 

In conclusion, we have reported a careful analysis of 
three neutron Bragg reflections that strongly constrain 
the character of the ordered phase in Tb2Ti2C>7 un- 
der the application of magnetic fields along [110]. The 
dominant scattering originates from Q ^ correlations, 
which are robust in fields as large as 9T, in contrast 
to the recent report of Cao et a/4 We have developed a 
general model based on the assumption that this scatter- 
ing has magnetic origin, and find that the AFM ordered 
moment is appreciable - conceivably on the order of 4.9 
\iB per spin if the system retains an Ising-like character. 
This represents a roughly threefold increase in the size 
of the AFM ordered moment relative to that refined in 
Ref H|. We contend that this discrepancy arises due to 
an extreme sensitivity to magnetic field misalignment. A 
possible test for this conjecture would be vector-magnet 
susceptibility measurements of the type reported in Rcf 
[23l We hope that this work motivates such further ex- 
perimental investigation, as well as theoretical work fo- 
cusing on the nature of the perturbation- induced ordered 
phases of this enigmatic magnet. 
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